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initio Quantum Dynamics and ab initio Molecular Dynamics 
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We report calculations of the dissociative adsorption of H2 
at Pd (100) covered with 1/4 monolayer of sulfur using quan- 
tum dynamics as well as molecular dynamics and taking all 
six degrees of freedom of the two H atoms fully into account. 
The ab initio potential-energy surface (PES) is found to be 
very strongly corrugated. In particular we discuss the influ- 
ence of tunneling, zero-point vibrations, localization of the 
nuclei's wave function when narrow valleys of the PES are 
passed, steering of the approaching H2 molecules towards low 
energy barrier configurations, and the time scales of the center 
of mass motion and the other degrees of freedom. Several "es- 
tablished" concepts, which were derived from low-dimensional 
dynamical studies, are shown to be not valid. 
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The presence of an adsorbate on a surface can pro- 
foundly change the surface reactivity. An understanding 
of the underlying mechanisms is of relevance for, e.g., cor- 
rosion, lubrication, and catalysis. Typically the breaking 
of molecular bonds is the rate-limiting step in a surface 
chemical reaction, and the model system is the dissocia- 
tion and adsorption of hydrogen. Whereas on clean Pd 
(100) the dissociation of H2 happens very efficiently, it 
is well known that a small amount of adsorbed sulfur re- 
duces (poisons) the surface reactivity significantly 
Density functional theory (DFT) calculations have shown 
that hydrogen dissociation on sulfur-covered Pd(100) is 
still exothermic, however, the dissociation is hindered by 
the formation of energy barriers in the entrance channel 
of the potential-energy surface (PES) 

We have recently extended this study and de- 
termined the six-dimensional PES of the system 
H 2 /S(2x2)/Pd(100) in great detail [j§ using DFT to- 
gether with the generalized gradient approximation 
(GGA) S . On an analytical representation of this ab ini- 
tio PES we have now performed six-dimensional quantum 
and classical dynamics calculations in which all hydrogen 
degrees of freedom are treated dynamically. The meth- 
ods are described in detail in Ref. [|| . The study is based 
essentially on the following approximations: (i) Born- 
Oppenheimer approximation, (ii) the substrate is kept 
rigid (Hi) the exchange-correlation functional is treated 
in the GGA, (iv) the DFT-GGA calculations are per- 
formed in the supercell approach, (v) the ab initio PES 
is fitted to an analytical form, 

Items i) and ii) mean that energy dissipation is not 



treated explicitly. This is appropriate because, while en- 
ergy dissipation plays a crucial role by taking away the 
adsorption energy, for the dissociation of H2 energy dis- 
sipation into substrate or electronic degrees of freedom 
is not significant. The reason lies in the involved time 
scales and the mass differences between hydrogen and S 
and Pd substrate atoms. 

The GGA treatment of the exchange-correlation func- 
tional [item Hi)] is the best treatment known to date. It 
leaves an uncertainty to total-energy differences typically 
below 0.1 eV (see, e.g., Ref. @). The error due to ap- 
proximation iv) was checked to be smaller than 0.1 eV 
|,|. The fit of the PES [item v)} by an analytic func- 
tion has been performed in such a way that the relative 
error between the fit and the original data, |Ai?|/|i?|, is 
smaller than 5 %. For further details of the DFT calcu- 
lations and the analytical representation of the PES see 
Ref. HJ. We note that the PES of this system is very 
corrugated, which represents a significant challenge for a 
proper treatment of the quantum dynamics. Even earlier 
empirical studies had not considered such corrugations 
(energy barriers between 0.09 eV for the optimum path- 
way and 2.5 eV for an approach over the S adatom pj). 
Nevertheless, we ensured that the error due to the finite 
number of channels considered in the quantum dynamical 
calculations is below 2 % of the sticking probability. 

The corrugation of the PES has several unexpected 
consequences. For example, it leads to zero-point vibra- 
tions (ZPV) in the quantum dynamics which reduce the 
quantum sticking probabilities compared to the classical 
results. Since, however, the ZPV energies in the single 
frustrated modes are rather small, the quantum particles 
do not realize completely the ZPV in these slow modes. 
Furthermore we find that steering of molecules to low- 
barrier configurations is effective at unusually high ki- 
netic energies. And in spite of the fact that the PES 
exhibits a early minimum barrier for dissociation, initial 
vibrational excitation of impinging molecules enhances 
the dissociation probability. These results confirm that 
results derived from statical analyses of the PES, as e.g. 
the barrier distribution, are not sufficient to obtain reac- 
tion probabilities, but a proper high- dimensional dynam- 
ical treatment of the dissociation process is essential. 

Figure |l| shows two cuts through the six-dimensional 
PES with the optimum pathway displayed in Fig. |l|a. The 
minimum energy barrier to dissociation (E\, — 0.09 eV) 
lies at a configuration where the bond-length of the 
molecule is is nearly unchanged compared to the gas 
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FIG. 1. Lines of constant energy for two cuts through the 
six-dimensional PES of H 2 at S(2x2)/Pd (100). The insets 
show the orientation of the molecular axis and the lateral Eh 
center-of-mass coordinates, i.e. the coordinates X, Y, 6, and 
4>. The coordinates in the figure are the H2 center-of-mass 
distance from the surface Z and the H-H interatomic distance 
r. Energies are in eV per H2 molecule. The contour spacing 
in Fig. |a is 0.1 eV, while in Fig. |l|b it is 0.2 eV. Fig. [i]a 
corresponds to the minimum energy pathway. 



phase value [|3|JJ]- In such a situation it is usually antic- 
ipated that the vibrational and translational degrees of 
freedom are almost uncoupled so that vibrational energy 
of the impinging molecules cannot be used to overcome 
the barrier (see, e.g., Ref. [HH). Consequently, it has 
been predicted [[| that the sticking probability of H2 at 
S(2x2)/Pd(100) should show no strong dependence on 
the initial vibrational state of the molecule. This predic- 
tion corresponds to the so-called Polanyi rules which have 
been formulated for gas-phase dynamics thirty years ago 
pof . We will show that these rules have to be modified 
for strongly corrugated PESs. 

The large repulsion for hydrogen dissociation closer to 
the sulfur atoms is illustrated in Fig. [l]b. While the disso- 
ciation path over the on-top position on the clean surface 
is hindered by a barrier of height 0.15 eV |11|, the ad- 
sorbed sulfur leads to an increase in this barrier height 
to 1.3 eV. 

Figure ^ compares our results for the sticking proba- 
bility as a function of the kinetic energy of the incident 
H2 beam with the experiment Q. In addition, also the 
integrated barrier distribution Pb(E) is plotted, which is 
the fraction of the configuration space for which the bar- 
rier towards dissociation is less than E and which is also 
called the "hole model" 

The calculated sticking probabilities are significantly 
larger than the experimental results. The agreement be- 
tween the adsorption experiment and our calculations 
is limited to the single aspect that adsorbed sulfur is 
a strong poison for the H 2 dissociation, and the onset 
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FIG. 2. Sticking probability versus kinetic energy for a H2 
beam under normal incidence on a S(2x2)/Pd(100) surface. 
Full dots: experiment (from ref. |lj); dashed-dotted line: inte- 
grated barrier distribution, which corresponds to the sticking 
probability in the hole model solid line: quantum me- 
chanical results for molecules initially in the rotational and vi- 
brational ground-state (QD); dashed line: classical molecular 
dynamics results for initially non-rotating and non-vibrating 
molecules (MD). The inset shows the quantum and classical 
results at low energies. 



energy for dissociative adsorption, calculated as E\ w 
0.12 eV, is in fact in excellent agreement with the ex- 
perimentally measured mean kinetic energy of hydrogen 
molecules desorbing from sulfur covered Pd(100) [ p^[ . 
Otherwise, in the energy range E x > 0.12 eV, the calcu- 
lated sticking probabilities are significantly larger than 
the experimental results. We consider this a most im- 
portant result, because progress in science is typically 
achieved by disagreement between theoretical and exper- 
imental results rather than by agreement. What are the 
likely reasons for the accounted disagreement? 

i) From the approximations assumed in the calcula- 
tions (see above) it is mainly the GGA which bears the 
uncertainty to be good enough for a reliable description 
of the topology of the PES. It should be noted that S(E{) 
probes the details of the energy corrugation of the PES 
and not just the lowest energy barrier. We note, however, 
that our results are not very sensitive to an uncertainty 
of the PES of < 0.2 eV, which we tested by adding some 
artificial functions to our ab initio PES, and that there 
is no indication in the literature (so far) that the GGA 
would not be able to deal with such a required, modest 
accuracy. 

ii) The other approximations have been carefully 
tested and we do not believe that they are relevant. 

Hi) Concerning the experimental results we note that 
the determination of the coverage is not without prob- 
lems. In fact, it is well possible that the adsorbate cover- 
age in the measurements was higher than reported. The 
sulfur adlayer was obtained by simply heating the sam- 
ple which leads to the segregation of bulk sulfur at the 
surface. The coverage of 8s = 0.25 was then identified 
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mode 


ZPV (eV) 


H-H vibration 


0.253 


polar rotation 


0.016 


azimuthal rotation 


0.013 


translation perpendicular to molecular axis 


0.027 


translation parallel to molecular axis 


0.027 


sum 


0.336 



TABLE I. Zero-point vibrations (ZPV) energies of the Eh 
molecule at the minimum barrier position. The H2 configu- 
ration corresponds to the situation of Fig. [j]a. The gas-phase 
ZPV energy of H 2 is \%u %3S = 0.258 eV. 



by monitoring the S132 and Pd33o Auger signals with re- 
spect to those at saturation coverage which is believed to 
be 6s = 0.5, and by the presence of a (2 x 2) LEED pat- 
tern. We speculate that the actual coverage at which the 
data of Fig. || were taken could well have been higher, 
e.g. due to some significant amount of subsurface sul- 
fur. Subsurface adsorbates can affect the surface reac- 
tivity significantly. We hope that more experiments of 
the type done by Rendulic et al. |l| will be performed 
and that also other techniques for producing the S ad- 
layer will be considered, in order to determine in what 
form S is actually present at the surface under different 
experimental conditions. 

We now turn back to the theoretical results. The clas- 
sical treatment of hydrogen dynamics over-estimate the 
sticking probability of H2 at S(2x2)/Pd(100) compared 
to the quantum results. For energies smaller than the 
minimum barrier height E\, the classical sticking proba- 
bility is of course zero, whereas the quantum results still 
show some dissociation due to tunneling, as the inset of 
Fig. U reveals. But for energies larger than E\> the clas- 
sical sticking probability rises to values which are almost 
50% larger than the quantum sticking probabilities. We 
have recently shown that ZPV can lead to a suppression 
of the quantum sticking probability compared to classi- 
cal results pjL4||. In Tab. | we have therefore collected 
the ZPV energies of the hydrogen molecule at the mini- 
mum barrier position. Since here the hydrogen molecular 
bond is still almost intact, the zero-point energy in the 
H-H vibration has only decreased by 5 meV from the gas- 
phase ZPV energy of H2. However, along the minimum 
energy path the molecule has to pass through a narrow 
valley of the PES in the lateral coordinates parallel to 
the surface due to the strong corrugation of the PES. 
This causes substantial ZPV energies in the frustrated 
lateral modes of the molecular center-of-mass of motion 
already relatively far away from the surface. Note that 
the ZPV energies in the frustrated rotational modes are 
fairly small in comparison. This is again due to the fact 
that the molecular bond is essentially not elongated at 
the minimum barrier position so that the molecular in- 
teraction with the surface is still rather isotropic. 

Because the sum of all ZPV energies at the minimum 



barrier position is 0.08 eV larger than the gas-phase ZPV 
energy of hydrogen, the effective minimum barrier height 
for the quantum particles should be increased by 0.08 eV 
to E^f — 0.17 eV. And indeed the quantum sticking 
probabilities are smaller than the classical sticking prob- 
abilities for energies larger than the classical minimum 
barrier height. The absence of ZPV in the classical dy- 
namics leads to the enhanced sticking probability in the 
classical calculations H]. However, the quantum stick- 
ing probability rises significantly already for energies well 
below the effective barrier height of 0.17 eV. This can be 
understood by considering the involved time-scales. A 
vibration with Jko = 50 meV, e.g., corresponds to a vi- 
brational period of 80 fs. The time it takes for an H2 
molecule with E\ — 0.12 eV to pass the barrier (width 
0.5 A) is only 25 fs. This means that during the cross- 
ing of the minimum barrier region the hydrogen molecule 
would hardly perform one complete vibration in any of 
the four frustrated modes given in Tab. Q. The concept 
that in the absence of tunneling a quantum particle needs 
at least the barrier energy plus the ZPV energies in order 
to cross a barrier is evidently not applicable for such slow 
modes. The quantum particles do not completely follow 
the zero-point motion perpendicular to the reaction path 
and thus do not additionally need the total ZPV energies 
in order to cross the barrier. 

The question arises why the difference between clas- 
sical and quantum sticking probabilities even increases 
with increasing energy although ZPVs should become less 
important at higher kinetic energies. An analysis of the 
ZPV shows that the sum of the four ZPV energies of the 
frustrated lateral and rotational modes actually rises to 
values of 0.2 eV along the reaction pathways. Although 
the H-H vibrational energy further decreases, the sum of 
all five ZPV energies becomes even larger than the value 
of 0.336 eV at the minimum barrier position. This sum 
plays a crucial role in the quantum dynamics not only at 
the minimum barrier position, but also later along the 
dissociation pathways. Since at different kinetic energies 
different parts of the barrier region with varying ZPV en- 
ergies are sampled, classical and quantum results are not 
just shifted by a fixed amount, and the ZPV need to be 
included in the full dynamical treatmemt. 

A comparison of the dynamical results with the inte- 
grated barrier distribution in Fig. || shows that the stick- 
ing probabilities are much larger than what one would 
expect from the hole model | fl2f . Our analysis of swarms 
of trajectories with different initial conditions shows that 
molecules directed towards high barriers are very effi- 
ciently steered to configurations and sites with lower 
barriers |l5|,|l6|]. Indeed the strongly repulsive potential 
above the sulfur atoms extends rather far into the gas- 
phase Thus the forces re-orienting and re-directing 
the molecules act on the molecules for a long time; there- 
fore they are so effective up to kinetic energies of more 
than 0.5 eV. At the clean Pd surface, on the contrary, 
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FIG. 3. Dependence of the quantum sticking probability 
versus kinetic energy for a hydrogen beam under normal inci- 
dence at a S(2x2)/Pd(100) surface on the initial vibrational 
state v. Solid line: Vi = 0, dash-dotted line: Vi = 1. The 
molecules are initially in the rotational ground state. 



for hydrogen dissociation on Pd(100), but large quanti- 
tive differences to the experiment exist. The huge corru- 
gation and anisotropy of the PES lead to ZPV in the 
quantum dynamics which reduce the quantum results 
compared to the classical ones. However, at the minu- 
mum barrier the additional ZPV correspond to very slow 
modes so that the quantum particles do not completely 
follow these modes. In addition, the huge corrugation 
and anisotropy cause large steering effects at unusual 
high kinetic energies. The multi-dimensionality of the 
PES furthermore leads to an enhancement of the stick- 
ing probability for vibrating molecules although the PES 
exhibits an early barrier for dissociative adsorption. 



at Ei sa 0.2 eV steering is no longer operative in the dis- 
sociation process |Tq| . This causes the surprising result 
that at E\ — 0.2 eV the sticking probability at the (2x2) 
sulfur-covered Pd(100) surface is almost as large as at 
the clean Pd(100) surface in spite of the fact that at the 
clean surface the fraction of open dissociation pathways 
is more than five times larger at this energy |L7| . 

The effect of initial vibrational motion on the dissoci- 
ation process is illustrated in Fig. || where state-specific 
sticking probabilities as a function of the incident kinetic 
energy are plotted. As noted above, since the PES has 
an early minimum barrier for dissociation, no coupling of 
the initial vibration to the dissociation is expected. How- 
ever, as Fig. |3| demonstrates, initial vibrational excitation 
leads to a significant increase in the sticking probability. 
This shows that the dissociation process cannot be un- 
derstood by simply analyzing the minimum energy path. 
As Fig. 0b reveals, there are molecular configurations for 
which the bond is significantly extended at the barrier po- 
sition. At higher kinetic energies the adsorbing molecules 
also probe such pathways along which the vibrational en- 
ergy can be efficiently used to overcome the dissociation 
barrier. The unexpected result of the vibrationally en- 
hanced dissociation in spite of an early minimum barrier 
to dissociation is actually in agreement with the experi- 
mentally observed vibrational over-population in thermal 
hydrogen desorption from sulfur covered Pd(100) Jl8| in- 
voking the principle of microscopic reversibility. We have 
also studied the dependence of the sticking probability on 
the initial rotational quantum number. The rotational ef- 
fects are similiar to the ones found at the clean surface 
and will be discussed in a forthcoming publication. 

In conclusion, we reported a six-dimensional dy- 
namical study of the dissociative adsorption of H2 at 
S(2x 2) /Pd(100) employing a PES obtained from detailed 
density functional theory calculations. The dynamical re- 
sults reproduce the poisoning effect of sulfur adsorption 
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